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Abstract. The Keck telescope's High Resolution Spectrograph (HIRES) has previously 
provided evidence for a smaller fine-structure constant, a, compared to the current labora- 
tory value, in a sample of 143 quasar absorption systems: Aa/a = (-0.57 ± 0.11) x 10" 5 . 
The analysis was based on a variety of metal-ion transitions which, if a varies, experience 
different relative velocity shifts. This result is yet to be robustly contradicted, or confirmed, 
by measurements on other telescopes and spectrographs; it remains crucial to do so. It is 
also important to consid er new possible ins trumental systematic effects which may explain 
the Keck/HIRES results. iGriest et alj J20CM arXiv:0904.4725vl) recently identified distor- 
tions in the echelle order wavelength scales of HIRES with typical amplitudes ±250 ms~'. 
Here we investigate the effect such distortions may have had on the Keck/HIRES vary- 
ing a results. Using a simple model of these intra-order distortions, we demonstrate that 
they cause a random effect on Aa/a from absorber to absorber because the systems are at 
different redshifts, placing the relevant absorption lines at different positions in different 
echelle orders. The typical magnitude of the effect on Aa/a is ~0.4 x 10~ 5 for individual 
absorbers which, compared to the median error on Aa/a in the sample, ~1.9 x 10~ 5 , is 
relatively small. Consequently, the weighted mean value changes by less than 0.05 x 10~ 5 
if the corrections we calculate are applied. Unsurprisingly, with corrections this small, we 
do not find direct evidence that applying them is actually warranted. Nevertheless, we urge 
caution, particularly for analyses aiming to achieve high precision Aa/a measurements on 
individual systems or small samples, that a much more detailed understanding of such intra- 
order distortions and their dependence on observational parameters is important if they are 
to be avoided or modelled reliably. 

Key words. Instrumentation: spectrographs - Techniques: spectroscopic - Cosmology: ob- 
servations - Quasars: absorption lines - Line: profiles 



1. Introduction 

The Standard Model of particle physics is 
parametrized by several dimensionless 'funda- 
mental constants', such as coupling constants 
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and mass ratios, whose values are not predicted 
by the Model itself. Instead their values and, 
indeed, their constancy must be established ex- 
perimentally. If found to vary in time or space, 
understanding their dynamics may require a 
more fundamental theory, perhaps one unify- 
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ing the four known physical interactions. One 
such parameter whose constancy can be tested 
to high precision is the fine-structure con- 
stant, a = e 2 /hc, characterising the strength of 
electromagnetism. Earth-bound laboratory ex- 
periments, conducted over several-year time- 
scales, which use ultra-stable lasers to com- 
pare different atomic clocks based on differ- 
ent ato ms/ions (e.g. Cs, Hg + , Al + , Yb + , Sr, 
Dy; e.g. lPrestage et al.ll 19951) . have limited the 
time-deriv ative of a to a/a = (— 1.6 + 2.3) x 
lO^yr 1 dRosenband et alj|2008l) 

Important probes of variations over cosmo- 
logical space- and time-scales are narrow ab- 
sorption lines imprinted on the spectra of dis- 
tant, background quasars by gas clouds asso- 
ciated with foreground galaxies (Bahcall et al. 
Il967l) . In particular, electronic resonance tran- 
sitions from the ground states of metallic atoms 
and ions are useful indicators of a variation for 
redshifts up to ~4 - when the universe was 
~10% of its current age - where the transi- 
tions are easily accessed from ground-based 
optical t elescopes. The many-multiplet (MM) 
method (iDzubaet all 1 199% I Webb et al.lll999h 
utilises the relative wavelength shifts expected 
from different transitions from different multi- 
plets of various atom/ions to measure a from 
quasar absorption spectra. The velocity shift, 
Av,, of transition i due to a small relative vari- 
ation in a, Aa/a <K 1, is determined by the 
^-coefficient for that transition, 



oj zJ = o) ,i + qt [(a z /a ) 2 - l] or 



Av; ^ ^ Aa qt 
c a o»oi 



(1) 



(2) 



where ojqj & a> z j are the rest-frequencies in the 
lab and at redshift z, ffo is the lab value of a and 
a z is the shifted value measured from an ab- 
sorber at z. The MM method is the comparison 
of measured velocity shifts from several transi- 
tions (with different ^-coefficients) to compute 
the best-fit Aa/a. FigureQ~]illustrates the wave- 
length shifts experienced by the transitions typ- 
ically utilized in MM analyses. 

Some evidence for a-variation has 
emerged over the last decade fr om quasar 
spectra observed with HIRES dVogt et al.1 



11994 on the Keck I 10-m telescop e in Hawaii. 
The fi rst tentative evidence in Webb et al. 
(1999) was subsequently strength ened with 
larger samples of absorbe rs dMurphy et al] 
2001a; iMurphv et all l2003l hereafter ImoUT 
MM analysis of 143 absorption spectra, all 
from the Keck/HIRES instrument, currently 
indicate a smaller a in the clouds at the 
fractional level Aa/a = (-0.57 + 0.11) x 10~ 5 

over the redshift ra nge 0.2 < Zabs < 4.2 

(IMurphv et alJl2004 hereafter lMOU) . 

Obviously, confirmation of such a po- 
tentially fundamental result must be made 
with many other telescopes and spectrographs. 
Similar MM studies using the Ultraviolet and 
Visual Echelle Spectrograph on the ESO Very 
Large Telescope in Chile are also beginning 
to yield constraints on Aa/a, but none yet 
rule out (or confirm ) the Keck/HIRES results 
dMurphy et al.ll20 08). At the same time, further 
searches for subtle systematic errors which, 
despite extensive searches (e.g. Murp hvet alj 
2001b), have evaded detection so far, must be 
considered in detail. Of particular importance 
are systematic errors in the wavelength cali- 
bration of the quasar spectra, which is usu- 
ally established using exposures of a stan- 
dard thorium-argon (ThAr) hollow-cathode 
emission-line lamp taken immediately before 
and/or after the quasar exposure. After the 
wavelength-pixel mapping is derived from the 
ThAr exposure, that solution is simply applied 
to the quasar exposure. Since the quasar and 
ThAr light illuminate the spectrograph slit dif- 
ferently and, in general, take slightly different 
paths through the spectrograph to the CCD, 
wavelength calibration errors may ensue. 

iGriest et alJ d2009h (hereafter lG09h re- 
cently observed a quasar absorption system us- 
ing Keck/HIRES with an iodine gas absorp- 
tion cell for more direct wavelength calibra- 
tion, thereby allowing comparison with the 
ThAr wavelength scale. For varying-a stu dies, 
the most important and robust result from lG09l 
is the identification of distortions in the ThAr 
wavelength scale across each echelle order 
(hereafter 'intra-order distortions'). See their 
figures 4 and 5. These may result from differ- 
ential and variable vignetting of quasar light 
compared to ThAr light, where only the former 
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Fig. 1. Transitions used in MM analyses versus their rest-frame wavelength and their sensitivity to varia- 
tions in a. Note the very different signature of a varying a in absorption systems where only the Mg and 
Fen transitions redwards of 2300 A are detected and fitted, compared with the more complicated signa- 
ture for systems containing a subset of the bluer transitions. Note that the relative strengths of the different 
transitions typically observed in quasar absorption systems is not portrayed here. 



enters HIRES from the telescope, the optical 
axes of which differ slightly, causing the quasar 
light p ath to rotate as the telescope tr acks the 
quasar dSuzuki et alJ2003tlG09h . lG0fl also dis- 
cuss velocity offsets between the ThAr and I2 
wavelength scales. If those results are robust, 
they are less important for varying a studies 
because, as equation (|2]i makes clear, such ve- 
locity offsets will not directly influence a Aa/a 
measurement 

This paper aims to assess the impact these 
intra-order distortions may have had on the 
Keck/HIRES results for varying a. The fol- 
lowing section describes the general effect the 
distortions will have and provides a crude cal- 

1 This is not strictly true when many quasar ex- 
posures are combined, as is generally the case for 
most of the Keck/HIRES sample. If different ve- 
locity shifts apply to the different exposures of the 
same quasar, and if the relative weights of the ex- 
posures vary with wavelength when forming the fi- 
nal, combined spectrum, then small relative velocity 
shifts will be measured between transitions at differ- 
ent wavelengths. Nevertheless, the effect on Aa/a of 
overall velocity shifts is of secondary importance to 
the intra-order distortions considered in detail here. 



culation of their expected magnitude. Section 
[3] details a more refined calculation of the 
correction to Aa/a for each absorber in the 
Keck/HIRES sample and Section [4] discusses 
the results. In section[5]we search for direct ev- 
idence for the need to apply the corrections and 
consider the effect of model errors on our cal- 
culations. We conclude in Section|6] 



2. The effect of intra-order 
wavelength distortions 

The reason intra-order distortions are problem- 
atic for varying a analyses is that they will 
produce velocity shifts between transitions. For 
an individual absorption system the transitions 
will, in general, fall at different positions along 
different echelle orders and will experience dif- 
ferent velocity shifts due to the wavelength cal- 
ibration distortions. The lG09l distortions gen- 
erally cause blueward shifts which increase to- 
wards the echelle order edges in either direc- 
tion from the center. This pattern seems to 
be generally repeated on each echelle order. 
However, because the I2 cell imprints absorp- 
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tion lines on the quasar spectrum only over 
the wavelength range ~5000-6200 A, no infor- 
mation about echelle orders outside this range 
can be obtained. Nevertheless, it is clear that 
for two absorbers at different redshifts, possi- 
bly with different transitions detected and fit- 
ted, different spurious shifts in Aa/a will be 
caused by the intra-order distortions. Indeed, in 
general, the effect on Aa/a will be random in 
sign and magnitude from absorber to absorber. 

A simple illustrative estimate of the ex- 
pected magnitude of the effect on Aa/a in a 
typical Mg/Fen absorber can be calculated as 
follows. Consider two transitions which, de- 
pending on the redshift of the absorber, fall 
at different positions along their echelle or- 
der. If w e mo del the intra-order distortions 
found by |G09| as a simple saw-tooth pattern 
of peak-to-peak amplitude Av amp = 500 ms -1 
then the root-mean-square (RMS) velocity dif- 
ference between the two transitions, averaged 
over redshift, will be A v amp / (2 y2) - 177ms -1 
if the transitions have arbitrary wavelengths. 
If one transition is from Mg and the other 
is from Fen, with typical rest wavelengths 
-2700A (i.e. frequencies -37000cm" 1 ), then 
the typical difference in ^-coefficients will be 
~ 1 250 cirT 1 . The resulting effect on Aa/a then 
follows from equation (0: |Aa/a colT | ~ 0.87 x 
10~ 5 . Adding more transitions will reduce this 
effect approximately as yN. Thus, with typi- 
cally 4-6 transitions observed in a single ab- 
sorber, the spurious effect the intra-order dis- 
tortions will have on its value of Aa/a will be 
|Aor/a colT | ~ 0.4 x 10~ 5 . 

Degradation in the accuracy of the wave- 
length solution near the order ed ges was al- 
ready considered a possibility in Web b et all 
(1999) . To test the possible effect of this 
IWebb et all artificially increased the 1-cr er- 
ror bar on Aa/a for absorption systems with 
transitions falling near order edges, thereby 
down-weighting those systems in the calcula- 
tion of the weighted mean Aa/a for the entire 
sample. The effect was found to be insignif- 
icant. However, a more detailed calculation is 
clearly warranted given the additional informa- 
tion about the particular form of the i ntra-o rder 
wavelength distortions identified by IG09I : we 
carry this out in the following section. 



3. Calculating the effect on the 
Keck/HIRES results 

In principle, the best way to understand 
the effect of intra-order distortions on the 
Keck/HIRES results would be to correct the 
wavelength scale of each extracted quasar ex- 
posure individually and then combine together 
exposures of the same quasar to again form 
the 1-dimensional spectrum. The^ 2 minimiza- 
tion of the Voigt profile fits to the absorbers 
could be run again on the new versions of the 
combined spectra, producing new, corrected 
values of Aa/a. This, clearly, would involve 
enormous effort and the results would still 
be subject to possible errors in our model of 
the intra-order distortions (hereafter 'model er- 
rors'). We therefore take a cruder approach, but 
one which still allows us to check our expecta- 
tions that (a) the correction to Aa/a should be 
random in sign and magnitude from absorber 
to absorber and (b) that the typical correction 
will be of order Aa/a cou « 0.4 x 10~ 5 as cal- 
culated above (Section|2]i. 

We first assume an appropriate functional 
form to counter the intra-order distortions of 
lG09t a velocity shift increasing from zero at 
the echelle order centre to +500 ms~' (i.e. a 
redward shift) at the order edges linearly in 
wavelength space. We consider a more com- 
plex model in Section [5] to illustrate the likely 
size of model errors in our results. 

For each absorption system in the 
Keck/HIRES sample, we must determine 
where each transition falls with respect to 
the order edges to estimate the velocity 
shifts between all the transitions due to the 
intra-order distortions. This is complicated 
by the fact that most quasar spectra used 
in the analysis were formed by combining 
extracted spectra of many separate exposures. 
For each absorption system, we derived the 
order positions for all relevant transitions from 
the extracted, wavelength calibrated spectrum 
for each quasar exposure. Thus, for each 
transition, i, an average velocity shift, Av;, was 
determined assuming that all quasar exposures 
contributing spectra of that transition did so 
with equal weight. While the signal-to-noise 
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ratio (SNR) did vary between exposures, this 
is a good approximation for most absorbers. 

From the velocity shifts between the tran- 
sitions in a given absorber, equation (f2]i can 
be used to compute the corresponding correc- 
tion to the value of Aa/a. That is, Aa/a corr 
is the slope (up to a factor of -2) of a lin- 
ear least-squares fit of the values of Av;/c ver- 
sus qi/ojoj. We use values for the labor atory 
frequencies, a>o, and ^-coefficients from |M03| 
in these calculations. However, two complica- 
tions arise here. The first is that not all tran- 
sitions have spectra of the same SNR. This is 
easily remedied by weighting the least squares 
fit with the squares of the SNR values mea- 
sured from the combined spectra in the con- 
tinuum around each transition. 

The second complication is that the shape 
of the absorption profile of each transition, and 
the relative strength of its constituent velocity 
components, also determine its relative con- 
tribution to constraining Aa/a in a given ab- 
sorber. For example, in an absorber with many 
strong, marginally resolved velocity compo- 
nents, a strong transition may appear very satu- 
rated. Thus, the centroids of the velocity com- 
ponents are very weakly constrained, if at all. 
But if several weaker transitions were also fit- 
ted, those velocity components will be opti- 
cally thin and their centroids will be strongly 
constrained in the weak transitions. Thus, in 
this example, the strong transition may not 
contribute much to the final constraint on Aa/a 
in this absorber. Properly taking this into ac- 
count is as 'simple' as re-running the^ 2 mini- 
mization of the Voigt profile fits to the absorp- 
tion profiles after correcting the spectrum of 
each transition i by its corresponding veloc- 
ity shift, Ay;. However, we have ignored this 
complication, thereby allowing a simpler, more 
illustrative calculation. Nevertheless, we have 
tested the importance of this simplification by 
re-running the x 2 minimization on several ab- 
sorption systems, at both low- and high-z, and 
find it to be fairly unimportant in most cases. 

4. Results 

Our estimates of the corrections to Aa/a 
for each of the 143 absorbers in the M04 



Keck/HIRES sample, calculated using the 
method described above, are given in Table [TJ 
These corrections should be added to the orig- 
inal values of Aa/a, which are also given in 
the table for convenience. The upper panel of 
Fig. [2] shows the corrections plotted versus the 
redshifts of the absorbers. Immediately we no- 
tice that the sign and magnitude of the cor- 
rections vary randomly from absorber to ab- 
sorber. Indeed, the mean correction, (-0.05 + 
0.05) x 10~ 5 , is consistent with zero, indicating 
that there is no strong influence on the aver- 
age value of Aa/a, as expected. The median 
magnitude of the corrections is 0.37 x 10~ 5 , in 
line with expectations from the very crude es- 
timate made in Section [2] And while the as- 
sumptions required to perform the calculation 
for each absorber, detailed in the last section, 
are not unimportant, they should not signifi- 
cantly affect these two main conclusions. 

It is important to emphasise that the sub- 
samples of absorbers below and above z a bs = 
1.8 are qualitatively different. All the absorbers 
at Zabs < 1.8 ('low-z') depend only on Mg and 
Fen transitions, whereas those at Zabs > 1-8 
('high-z') tend not to include Mg and generally 
contain the wider variety of transitions, with a 
mixture of ^-coefficients, bluewards of A mst = 
2300 A in Fig. [TJ As illustrated in Fig.[Q the 
signature of a varying a for Mg/Fe n absorbers 
is therefore very simple, while that for the 
higher-z systems is more complicated, strongly 
depending on which transitions are detected 
and fitted. In this sense, absorbers with only 
Mg/Fe ii transitions fitted are more susceptible 
to simple instrumental systematic errors which 
cause long-range, low-order distortions of the 
wavelength scale and also astrophysical effects 
which might shift velocity components of Fe n 
relative to those of Mg. However, for the low-z 
Mg/Fe ii absorbers in the Keck/HIRES sample, 
the scatter around the weighted mean Aa/a 
was c ons istent with the individual error bars in 
M03 and lM04l This remains true even after the 
corrections for intra-order distortions are ap- 
plied: x 2 = 76.03 around the weighted mean 
Aa/a for the 77 low-z absorbers. 

The same canno t be said of the high- 
z absorbers: in IM03I and|M04j we identified 
a sub-sample of high-z absorbers for which 
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Table 1. Values of Aa/a and corrections for possible intra-order distortions of the wavelength scale. The 
first three columns specify the B1950 quasar name, the emission and absorption redshifts. The fourth col- 
umn provides the transitions used to determine Aa/a in each absorber; each letter represents a transition 
and the key is provided in table 2 of M03. The fifth column is the raw value of Aa/a together with the 
formal 1-cr statistical error. The systematic error term, 6(Aa7tt) sys , is non-zero only for "high-contrast" 
absorbers; see text for explanation. The seventh and eighth columns provide our estimate of the correc- 
tion, Aa/a cmn for intra-order distortions, to be added to Aa/a, and the systematic error term, S(Aa/a) c s °l r , 
which was calculated in the same way as 6(Aa/a) sys but using the corrected values of Aa/a. The final 
column specifies the observational sub-samples employed: A = "Previou s low -z sample" from lM03l ; Bl 
- "P revious high-z sample" from lM03l : B2 = The 15 absorbers added in M04; C = "New sample" from 
M03. Only a excerpt of the table is provided here; the complete table is available from the authors or 
http://astronomy.swin.edu.au/~mmurphy/pub.html . 



B 1950 name 




Zabs 


Transitions 


Aa/a 


6(Aa/a) sys 


Aa/a con . 


6(Aa/a)^ 


Sample 










[io- 5 ] 


[IO" 5 ] 


[io- 5 ] 


[io- 5 ] 




1634+7037 


1.34 


0.99010 


bcnpqr 


1.156 ±2.399 


0.000 


0.498 


0.000 


A 


0019-1522 


4.53 


3.4388 


ghl 


0.937 ±3.912 


0.000 


1.414 


0.000 


Bl 


0100+1300 


2.68 


2.3095 


efgjklmvw 


-3.949 ± 1.370 


1.754 


0.063 


1.707 


Bl 



the scatter in Aa/a exceeded expectations 
based on the individual errors. In these 27 ab- 
sorbers both very strong and very weak transi- 
tions were fitted simultaneously - i.e. they are 
"high-contrast" absorbers - making the multi- 
component Voigt profile fitting process more 
difficult and error-prone. These systems could 
be 'under-fitted' - too few velocity compo- 
nents used to model each absorber - and this 
may lead to systematic errors in individual ab- 
sorbers which are random in sign and mag- 
nitude. We demonstrated the ef fect of under- 
fitting with simulated spectra in Murphv et al. 
(2008). Just as in lM03l and lM04 the aver^ 
age systematic error in the sample of 27 high- 
contrast absorbers is that which, when added 
in quadrature to the individual statistical errors, 
reduces^ 2 per degree of freedom,^- 2 , to unity 
around the weighted mean Aa/a for those ab- 
sorbers. The fact that we find a very similar 
systematic error term using the corrected val- 
ues of Aa/a, 1.71 x 10~ 5 , as found from the 
uncorrected values, 1.75 x 10~ 5 , indicates that 
the extra scatter in Aa/a for these absorbers 
does not arise from the intra-order distortions; 
our postulate that it arises from the mixture 
of strong and weak transitions fitted, and the 
resulting under-fitting, remains. These system- 
atic error components are shown for the high- 
contrast absorbers in Table [TJ 



The middle panel of Fig. [2] shows the cor- 
rected values of Aa/a with their 1-cr statistical 
errors and systematic error components added 
in quadrature. The lower panel shows a binned 
version of these results, where the weighted 
mean Aa/a of the 13 absorbers in each bin 
is shown with its 1-cr e rror. As with the un- 
corrected results in |M04 Aa/a is consistently 
smaller a in the absorption systems compared 
to the current laboratory value. Table|2]quanti- 
fies this: the weighted mean Aa/a for the full 
sample, Aa/a = (-0.61 ± 0.11) x 10~ 5 , differs 
only slightly from the uncorrected value from 
M04, (-0.57 ± 0.11) x 10~ 5 , as expected. 

Table|2]gives the statistics for th e low- and 
high-z sub-samples. It was shown in M03| that 
the low- and high-z samples respond, on av- 
erage, in opposite ways to simple, long-range 
distortions of the wavelength scale. Note that 
the weighted means for the low- and high-z 
absorbers are similar and both depart signif- 
icantly from zero, even after the corrections 
are applied. This also quantifies the fact that 
the evidence for a varying a is dominated by 
the low-z absorbers, with the evidence at high- 
Z being at the ~G-cr level both before and af- 
ter correction. Still, as mentioned above, if 
long-range wavelength calibration distortions 
remain in the data, the low- and high-z sam- 
ples' weighted mean Aa/a values should have 
opposite sign. This is an important internal 
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Fig. 2. Corrections to Aa/a and corrected values of Aa/a versus absorption redshift and cosmological 
look-back time (H = 71 km s~ Mpc~', fl m0 = 0.27, fl A0 = 0.73). Top: Estimated corrections to Aa/a due 
to intra-order distortions. Note that their sign and magnitude vary randomly from absorber to absorber. The 
median magnitude of the corrections is 0.37 x 10~ 5 . Middle: The corrected values of Aa/a with their error 
bars representing the quadrature sum of the 1-cr statistical and systematic errors. Bottom: Weighted mean 
values of Aa/a, with 1-cr errors, in bins containing 13 absorbers each. 



Table 2. Basic statistics for the entire sample and redshift sub-samples. The second column gives the 
number of absorbers, N.^, in each sample. The statistics before any corrections for intra-order distortions 
are applied are in columns 3, 4 & 5, while the remaining columns provide the same statistics after the Aa/a 
values are corrected. (Aa/a) w is the weighted mean and (Aa/a) is the unweighted mean. 







Before correction 




After correction 




Sample 


A'abs 


<Aa/or) w 


{Aa/a) 


Median 


(Aa/a) v 


{Aa/a) 


Median 






[10- 5 ] 


[io- 5 ] 


[io- 5 ] 


[io- 5 ] 


[io- 5 ] 


[io- 5 ] 


Fiducial 


143 


-0.57 ± 0.11 


-0.53 ± 0.19 


-0.42 


-0.61 ± 0.11 


-0.58 ± 0.20 


-0.56 


Zabs < 1-8 


77 


-0.54 ±0.12 


-0.61 ± 0.22 


-0.39 


-0.58 ±0.12 


-0.60 ± 0.23 


-0.58 


Zabs > 1-8 


66 


-0.74 ± 0.27 


-0.45 ± 0.34 


-0.43 


-0.79 ± 0.27 


-0.57 ± 0.34 


-0.55 



consistency check that is only available when 
one compares low-z Mg/Fe n with higher-z ab- 
sorbers constraining a greater diversity of tran- 
sitions. Of course, we must also recognise that 



absence of evidence for systematic errors ca- 
pable of explaining the consistency of the low- 
and high-z Aa/a values is not evidence for 
their absence. 
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5. Are the corrections warranted? 

We saw in the previous section that correcting 
the Aa/a values for in tra-order distortions of 
the kind found bv lG09l makes no difference to 
the overall conclusions, nor for some more de- 
tailed aspects of the Keck/HIRES results. So, 
having calculated reasonable estimates of the 
corrections, can we find evidence that applying 
them really is warranted? 

For example, if the corrections are impor- 
tant, we should expect the distribution of Aa/a 
values around the mean to significantly narrow 
after applying the corrections. We should also 
expect a significant anti-correlation between 
the values of Aa/a and the corrections. We find 
neither of these effects. 

Table [2] provides the mean Aa/a and its 
standard error - i.e. the RMS/ -yW a bs - for the 
full sample and low- and high-z sub-samples. 
Only very small changes in the RMS are ev- 
ident in each case, indicating that the correc- 
tions do not remove a significant amount of 
scatter in the Aa/a values. A small Monte 
Carlo simulation can be used to gauge how 
much difference in the RMS we should expect 
if the corrections were important. We gener- 
ated random absorber data-sets, with the same 
size, total errors and corrections as the real 
Keck/HIRES sample. The corrected Aa/a val- 
ues in each realisation were set to -0.6 x 
10~ 5 and then randomised, according to the 
Gaussian errors for individual absorbers. The 
RMS of this sample, and the same realisation 
with the corrections removed from the Aa/a 
values, were compared. Differences in RMS 
values between the corrected and uncorrected 
realizations of > 0.1 and > 0.2 x 10~ 5 occurred 
37 and 4 % of the time by chance alone. Thus, 
the lack of RMS differences between our real 
corrected and uncorrected samples cannot be 
used as evidence that the corrections are not 
meaningful. But the Monte Carlo simulation 
does indicate that the corrections are too small, 
in comparison to the total errors on Aa/a in in- 
dividual absorbers, to make a large difference 
to the sample overall. 

Similarly, the Spearman rank correlation 
and Kendall's t tests find insignificant anti- 
correlations between Aa/a and the corrections 



for the full sample or sub-samples. Nor are ab- 
sorbers with large Aa/a errors masking an un- 
derlying anti-correlation; using only absorbers 
with total errors (quadrature sum of statisti- 
cal and systematic errors) less than 3 x 10~ 5 , 
gives similarly insignificant results. However, 
again a small Monte Carlo simulation confirms 
that, effectively, this is not unexpected given 
the small size of the corrections relative to the 
larger total errors on individual Aa/a values. 

To summarise, while we do not find evi- 
dence that the corrections we calculate need 
to be applied, it is their relative smallness in 
general which precludes a clear test for this. In 
essence, the formal errors in the Keck/HIRES 
sample still dominate over the potential sys- 
tematic errors caused by the intra-order distor- 
tions identified by |G09l 

However, an important assumption so far 
has been the particular form of in tra-order dis- 
tortion we have used; the lG09l results have 
been approximated with a simple saw-tooth 
patte rn in every echelle order. Because the 
G09 analysis relies on 12-cell calibration, it can 
only probe the intra-order distortions over the 
wavelength range ~5000-6200A. Even over 
that relatively short wavelength range, there is 
some evidence for a slow decrease in the peak- 
to-peak amplitude of the saw-tooth pattern in 
bluer orders. How this extrapolates to longer 
and shorter wavelengths than 6200 and 5000 A 
(respectively) is unknown. Thus, model errors 
may well exist in our estimates of the correc- 
tions to the Keck/HIRES Aa/a values above. 

To test the importance of this, we can re- 
peat the calculation of the corrections using 
a somewhat different model of the intra-order 
distortions: we fix Av amp to 500 m s~'at 5500 A 
and increase (decrease) it linearly with slope 
0.55 ms~' A -1 above (below) 5500 A and en- 
force a maximum (minimum) amplitude of 
1000 ms -1 (200 m s~') at redder (bluer) wave- 
lengths outside the range covered by the I2 cell 
calibration. The RMS difference between the 
old and new corrections is 0.56 x 10~ 5 and, 
with the new corrections, the weighted mean 
becomes (-0.62 + 0.1 1) x 10~ 5 over the whole 
sample and (-0.61 ± 0.12) x 10~ 5 [(-0.66 + 
0.27) x 10~ 5 ] at low-z [high-z]. Comparison 
with the values in Table [2] reveals that the 
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model errors in this case are very small. Of 
course, this model is not very different to our 
original, simpler one. For example, it may be 
that the intra-order distortions have a com- 
pletely different shape, phase and/or amplitude 
for different quasar observations. This possibil- 
ity must be explored with future observations. 

6. Conclusions 

Exploring systematic effects which may ex- 
plain the Keck/HIRES evidence for a varying a 
is clearly an important problem. However, vari- 
ous properties of the Keck/HIRES results make 
the task difficult; the consistency between the 
average Aa/a values in the low-z Mg/Fe n and 
the more diverse high-z systems being an im- 
portant one. We have demonstrated here that 
if we model the intra-order distortions iden- 
tified by lG09l in a simple way, and extrapo- 
late the model to all echelle orders (not just 
those within the 12-cell calibration range) the 
effect on the overall evidence for varying a is 
very small. As expected, the distortions affect 
individual values of Aa/a randomly in sign 
and magnitude, because the differing redshifts 
place the transitions at varying positions with 
respect to echelle order edges where the dis- 
tortions are worst. Indeed, the effect for a typ- 
ical absorber is to shift Aa/a by 0.4 x 10~ 5 
which, compared to the median error on Aa/a, 
~1.9 x 10~ 5 , is too small for us even to find 
direct evidence of the need to apply the correc- 
tions we calculate. 

Despite the above conclusions, it is im- 
portant to emphasise that we do not yet fully 
understand the origin of t he intra-order dis- 
tortions identified by IG09L how they depend 
on various observational parameters (e.g. tele- 
scope pointing direction, temperature, time 
etc.) and therefore how they may differ- 
ently affect spectra of different quasars in the 
Keck/HIRES sample. Indeed, lG09l show that 
separate exposures taken through the I2 cell 
seem to have somewhat different intra-order 
distortion patterns. And while we have made 
a simple attempt to address this problem of 
model errors, it is not enough to completely 
rule out intra-order distortions as an important 
systematic error for the Keck/HIRES results. If 



future Keck/HIRES observations allow much 
smaller statistical errors on Aa/a in individual 
absorbers or small samples, i ntra-o rder distor- 
tions like those identified by lG09l must either 
be eliminated or carefully modelled. 
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